Slave-boson field fiuctuation approach to the extended Falicov-Kimball model: 
charge, orbital, and excitonic susceptibilities 
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Based on the S'0(2)-invariant slave-boson scheme, the static charge, orbital, and excitonic sus- 
ceptibilities in the extended Falicov-Kimball model are calculated. Analyzing the phase without 
long-range order we find instabilities towards charge order, orbital order, and the excitonic insulator 
(EI) phase. The instability towards the El is in agreement with the saddle-point phase diagram. We 
also evaluate the dynamic excitonic susceptibility, which allows the investigation of uncondensed ex- 
citons. We find qualitatively different features of the exciton dispersion at the semimetal-EI and at 
the semiconductor-EI transition supporting a crossover scenario between a BCS-type electron-hole 
condensation and a Bose-Einstein condensation of preformed bound electron-hole pairs. 

PACS numbers; 71.30.-|-h, 71.35-y, 71.35.Lk, 71.28.-|-d 
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I. INTRODUCTION 

At low temperatures electronic correlations can cause 
anomalies at the semimetal-semiconductor (SM-SC) 
transitionji Half a century ago, Motlj^ argued that in 
a SM with a very low carrier density the Coulomb at- 
traction between electrons and holes should lead to the 
spontaneous formation of electron-hole bound states (cx- 
citons), and the system would become insulating. Shortly 
afterwards, Knojs^ noticed that a SC is unstable against 
the spontaneous formation of excitons if the exciton bind- 
ing energy overcomes the gap energy separating valence 
and conduction band. Both arguments suggest a new 
distorted phase, an exciton condensate known as the ex- 
citonic insulator (EI), to be the crystal ground state. The 
SM-EI transition is mathematically similar to the BCS 
theory of superconductivity, while the SC-EI transition 
can be treated as a Bose-Einstein condensation (BEC) of 
preformed excitons. Hence, the EI is discussed in view of 
a BCS-BEC crossover scenario in a solid^"— 

Whilst theoretically predicted a long time ago^ (for re- 
cent reviews see Ref. |9|) no conclusive experimental proof 
of the existence of the EI has been achieved yet. How- 
ever, there are a few promising candidates. In the mixed 
valence compound TmSeo.45Teo.55 detailed studies of the 
pressure-induced SC-SM transition suggest that excitons 
are created in a large number and condense below 20 
KiiS More recently, several transition-metal dichalco- 
genides were reported to exhibit an EI phase. Angle- 
resolved photo emission spectra (ARPES) measurements 
of Ta2NiSe5 traced back the extreme valence band top 
flattening at low temperature to an EI ground stateiii 
ARPES data of ir-TiSc2 indicate that the EI is the driv- 
ing force for the charge-density-wavc (CDW) transition 
in this materialJ^ 

From a theoretical point of view, the description of the 
EI with a Falicov-Kimball-type model seems promising. 
The original Falicov-Kimball (FKM) modeli^ contains 
itinerant c-electrons (with bandccnter and hopping 




FIG. 1: (Color online) Hartree-Fock ground-state phase dia- 
gram of the EFKM in two dimensions for Coulomb strength 
(7 = 2. The difference between CDW and SOO is explained 
in the text. The black solid line represents the second-order 
transition from an EI to a BI (band insulator), the dashed 
line represents the first-order CDW/SOO-EI transition. The 
label 0"*" emphasizes that the EI phase is present only for finite 
/-bandwidths. 



amplitude t^) that interact via a local Coulomb repul- 
sion U with localized /-electrons (with energy level Ef)^ 
where the spin is neglected. Since the local /-electron 
number is strictly conserved in the FKM, /-c-coherence 
cannot be established One way to overcome this short- 
coming is to include an f-c hybridization^^ As shown in 
Refs. flaflTl . the extension by a finite /-bandwidth also 
induces f-c coherence. The model with a direct /-/ 
hopping (with hopping amplitude tf) is called the ex- 
tended Falicov-Kimball model (EFKM) and has previ- 
ously been used to describe different properties of the 
EI phase i^i^ d^d^ The ground-state phase diagram of the 
EFKM was determined with a constraint path Monte 
Carlo (CPMC) technique for one and two dimensions 



2 



(ID and 2D) in the strongi^ and intermediate coupling 
regimeil as well as in the Hartree-Fock (HF) approx- 
imation for 2D,— 3D,— and infinite dimensions.— All 
approaches yield a qualitatively similar phase diagram. 
Figure [T] displays the HF ground-state phase diagram for 
[/ = 2 in 2D, depicting the generic situation. It was 
shown previously that Fig. [T] agrees with the CPMC data 
even quantitativcly>21i Besides full c- and full /-band in- 
sulator (BI) regions, the EFKM ground-state phase di- 
agram exhibits three symmetry broken phases: the EI, 
a CDW, and a staggered orbital order (SOO). The EI 
is characterized by a nonvanishing average (c^/). The 
CDW is described by a periodic modulation in the total 
electron density comprising both /- and c-electrons. The 
SOO is characterized by a periodic modulation in the 
difference between the /-electron and the c-clectron den- 
sity, which may be accompanied by a CDW. The SOO 
(CDW) establishes the ground state for the symmetric 
case {Ef = Ec) for aU ratios of —tf/tc (for the CDW 
the point \tf\ = \tc\ has to be excluded, see below). 
These phases are rapidly suppressed in favor of the EI if 
Ef 7^ Ec- Between the uniform EI phase and the CDW 
or SOO phase there is a first-order phase transition. The 
EI-BI transition is of second order. For tf —Q the FKM 
is recovered, and the EI phase cannot be realized. 



II. THEORY 



A. Model Hamiltonian 



Expressing the orbital fiavor by a pseudospin variable 

a =t, i, where c||' = /f^^ and c^f = c|^\ the EFKM can 
be written as an asymmetric Hubbard model. 



H 



{Ea 



(1) 

where c^P annihilates (creates) a cr-band electron at the 
Wannier site i and n^^ — c\^c^^ is the corresponding 
number operator. E„ denotes the bandcenter of the a- 
electron band, ^ gives the chemical potential, t„ is the 
hopping amplitude, and U measures the Coulomb inter- 
action strength. In what follows we consider i?^ = 0, 
< 0, = 1, and < 0. All energies are mea- 
sured in units of t^. We restrict ourselves to tit^ < 0, 
i.e., the valence band top and the conduction band min- 
imum are located at the Brillouin zone center. Moreover 
we exclusively investigate the half-filled band case, i.e.. 



N 



1, where N is the number of lattice sites. 



B. Slave-boson functional integral representation 



For the investigation of electron correlation effects the 
Gutzwiller approximation^^ is an established technique. 
Kotliar and Ruckenstein introduced a scalar slave-boson 
(SB) scheme which reproduces the Gutzwiller solution 
of the Hubbard model as a saddle-point 3^ A manifestly 
spin-rotation invariant form of the SB representation has 
been worked out for the Hubbard model^i and for multi- 
band Hubbard modelsi^ We have developed an 5*0(2)- 
invariant SB approach for the EFKM (Ref. [l^ that re- 
produces the HF result for the EI phase boundary at 
T = 0, but leads to a substantial reduction of the crit- 
ical temperature. It is the aim of this work to include 
Gaussian fluctuations around the saddle-point^"— in the 
50(2)-invariant SB scheme at zero and finite tempera- 
ture. This offers an opportunity to calculate suscepti- 
bilities for investigating instabilities against long-range 
ordered phases and the formation of excitons. 



The paper is organized as follows. In Sec. [IT] the model 
Hamiltonian and the SB scheme are introduced. More- 
over the saddle-point approximation is given and the cal- 
culation of response functions within the SB scheme is 
explained. In Sec. IIIII we present numerical results for 
the instabilities toward the CDW, the SOO, and the EI 
phase. Finally we investigate the formation of excitons 
in the phase without long-range order. Section |IV] sum- 
marizes our results. 



Following Refs. Il9ll24l the Hilbert space is enlarged by 
introducing auxiliary bosons: e|^^ , related to an empty 

site, d[^\ related to a doubly occupied site, and p'^.^\ 
related to a singly occupied site. 



|0,) -> e^lvac), 
p 



(2) 
(3) 
(4) 



where |vac) means the vacuum state. The matrix opera- 
tor p^^' is given as 



^(t) 



1 



"if 

2 ^ 2 \ „(t) , j„(t) „(t) _ (T) 



1 /r,(t) , Jt) Jt) 
(t) _ 1 / PiQ + Piz Pix ~ V 



(5) 

where /i = 0,x,y, z. Tq denotes the unit matrix and 
T = (t^, Tj^, T^)"^ is the vector of the Pauli matrices. The 
fermionic degrees of freedom are captured by the pseud- 
ofermions c| = (cl^.cl) and = (5,.^, 5,-^)^. 

Unphysical states of the extended fermion-boson Fock 
space are excluded by two sets of local constraints. 



(1) 



ele, -K2TrpV + djci, -1=0, 



(6) 



(2) 



e,6l+2£tp^ + 4d,ro-ro = 0. (7) 



Since the bosonic occupation number of one site is 
coupled to the fermionic occupation, the bosons have to 
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change simultaneously when an electron is created or an- 
nihilated. This is achieved by introducing the bosonic 
hopping operator Zj, 



where t = ^ ^ . The trace in Eq. ([Ti|) extends over 

time, space, and spin variables. For the half-filled band 
case Eqs. ([T4|) and (|T5|) are an exact representation of the 



(8) partition function of the EFKM. One obtains 



The choice of Zj is not unique. We choosei2i2i 

Z, = LAM^P^K^ + L,gM,d,N^ 

with 

L, = [(l-4rf,)ro-2£tpJ-i/2, 
= [l + ele,-H4d„ + 2TrpV]i/2 



(9) 

(10) 

(11) 
(12) 



and p[^Jp, — PP'Pi-p'-p 1 which guarantees the correct 
free-fermion result on the mean-field level. The Hubbard 



interaction term is bosonized via 



The resulting coupled fermion-boson system is evalu- 
ated within a functional integral representation. Then, 
the bosons become complex fields and the fermions are 
given by complex Grassmann fields. The Lagrange mul- 
tipliers A-^\ A-g\ x[l\ x[y\ and A-f^ are introduced to 
enforce the constraints (O and ([7]). Exploiting the gauge 
symmetry of the action and transforming the Lagrange 
multipliers into real time-dependent Bose fields we can 
remove the phases of pto, piz, and e;. Using the Grass- 
mann integration formula, we obtain the grand canonical 
partition function given by a functional integral over Bose 
fields only. 



Z = J D[e]D[po]D[p*,pjD[p;,Py]D[p,]D[d\d] 

i?[A(i']D[A^'^]D[A(2)]e-^ (13) 
with the effective bosonic action 

[dr\Y.[-X^^+>^^ei + j:i^^^>^)\p.,\^ 
•i ^ i II 







r(2) 



Tr In 



-PioiPi 

+ PixdrPix + PiydrPiy 

\-G7.\ At,t')\, 



d*drd. 



(14) 



where Pi = ipix,Piy,Piz) and Xf ^ = (A|^\ A,[^^\ A^^f^). The 



inverse Green propagator is given by 
-"-1 / i\ _ ( a 1 .. \(2 



6,^ 5{r - r') 
(15) 



C. Saddle-point approximation 

To proceed we approximate all bosonic fields by their 
time-averaged values (static approximation), i.e., the 
bosonic fields are taken to be real. Moreover, we look 
for uniform solutions, that is, the Bose fields are taken 
to be independent of the lattice site. 

We restrict ourselves to the phase without long-range 
order, which we denote as paraphase. The saddle-point 



equations for the paraphase {px 
are 



Py 



A 



(2) 



A. 



(2) 



POPz 



\ + \ln^ni{l - z2) 



1 

2^ 



Al^ 



U 



-2po^z^ (2 

Pz f 
^Po \2d^ 
2rf2 ^2 



Pi 



Pi 
,2^2 



pI) + z'^pI+pI 



z2e(0) 



-^3±^zM0)-2A(2)^ 
pIcP Po 



N 

k 
k 



where 



rika = [exp{f3Ek^) + 1] ^ , 



P 



P 



A, 



(2) 




= 0) 

(16) 
(17) 

(18) 
(19) 
(20) 
(21) 

(22) 



(23) 
(24) 
(25) 



On a D-dimensional hypercubic lattice, 7k = 
2^j^-^cosfci. The chemical potential is determined by 
the condition 



k,CT 



(26) 



The quasiparticle gap Eg indicates the splitting of the 

t- and |-band (in the paraphase) , which is caused by the 

(2) 

correlation-induced quasiparticle bandshift Az . For a 
D-dimcnsional hypercubic lattice. Eg is given by 



En 



\E,\ + \2Xf^\-2DzWt^\ + \t,\) 



(27) 
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For a SM, Eg < and for a SC, Eg > 0. 

We obtain the EI phase boundary by solving the SB 
gap equation, 



1 = 



1 



1 Jikt - ni^l 



PoPz 



N^E, 



kt 



E^ 



(28) 



rcsuhing from Eqs. (63) and (65) of Ref. [T^ 

The gap equation (f28| captures both the BCS and the 
BEG situation, but it cannot discriminate between them. 
To this end, we follow an idea from Ihle et al. (Ref.Q) and 
investigate the excitonic susceptibility in the paraphase. 



D. Gaussian fluctuations 

In order to study response functions, we take into ac- 
count Gaussian fluctuations around the saddle point for 
the paraphase, that is, $i„ = + (5$,:„, where $,„ = 



{ ^2 ; PiQ 7 Pix 7 Piy i P 

the action is given by 



J \(2) x(2) ,,(2) 



X^J}. Then, 



(29) 



q,a,h 



where the bar denotes the saddle-point value. 

In order to achieve comparability with the saddle-point 
results, we start the fluctuation calculation from the same 
level of approximation as for the saddle-point calculation, 
i.e., we first perform the static approximation and con- 
sider only the fluctuations of the 11 real- valued fields $io- 

The fluctuation matrix can be calculated according to 



Sab{q,q') 

— y 

MR ^ 



-iqRi 



2Np 



Sab{q)5q-q' ■ 



-iq' Rj 



(30) 



Here, we use the shorthand notation Ri — (Ri,T) and 
q = (q, Wm), where r is the imaginary time, uim = iirm/ (i 
denote the bosonic Matsubara frequencies, R,; is the po- 
sition vector, and q is the wave vector. 

The response functions can be expressed in terms of 
the SB field fiuctuations using the local constraints (|6]) 
and ([7|). The charge susceptibility reads 

Xc{q) = {5[n^{-q) + ni{-q)]5[n^{q)+n^{q)]) 
= 4{e^{Se{-q)Se{q)) -2ed{Se{-q)Sd{q)) 

+d^{Sdi-q)Sdiq))) . (31) 

The orbital susceptibility is given by 

Xo{q) = {6[n^{-q) - n:^,{-q)]S[n^{q) - ni{q)]) 

= 4{pl{Spo{-q)Spo{q)) +2p,po{Spo{-q)Sp,{q)) 
+pl{Sp,{-q)dpM)) ■ (32) 



Gonsidering the creation operator of an onsite electron- 
hole pair— 



bl 



1 



(33) 



the electron-hole susceptibility, hereafter denoted as ex- 
citonic susceptibility, is given by 

XX (q) = {Sb^SbD 

= Po [{SPx{-q) Spx{q)) + {Spy{-q) Spy{q)) 
-~i{Spy{-q) 6px{q)) + i{6px{-q) Spy{q))] . 

(34) 

The correlation functions may be expressed as func- 
tional integrals over Bose fields: 

{s^a{-q)SMq)) = I y Dm 6^ahq)SMq) ^'^^"^ ■ 

(35) 

Hence, the correlation functions are related to the inverse 
fiuctuation matrix by 



1 



{s<Pai-q)SMq)) = ^sjiq) 



(36) 



It turns out that for the paraphase the 11 x 11 fluctu- 
ation matrix decomposes into a 7 x 7 matrix containing 
the charge fluctuations {Se, SpQ, Sd, S\^^\ '^^o'^'') ^-^^^ the 

(2) 

orbital fluctuations {Spz, SXi ') and into a 4 x 4 matrix 
containing the electron- hole pair fluctuations {Spx, SXx , 

Spy, (5 Ay ). The SB fields are obtained by solving the 
saddle-point equations ((T5)) - (P^ self-consistently. 

The description of the GDW and SOO requires the in- 
clusion of inhomogeneous solutions with a periodic mod- 
ulation in the densities, {ui^} = rio- -I- ^o- cos(QRj), where 
the order vector in 3D is given by Q = (7r,7r,7r). The 
GDW and SOO order parameters are i5cdw = ^(i^t + ^i) 
and 6soo = ^{S^ - ^i), respectively^ If \S^\ ^ 
the GDW and SOO describe the same symmetry broken 
state. We can investigate the formation of both phases 
without generalizing the SB formalism to a bipartite lat- 
tice by calculating the static (w = 0) charge and orbital 
susceptibility with order vector q = Q, given by 

xc - Xc(Q,o) 

= 2[e2(5-i)ee + rf'(5-i)dd-2ed(5-i)ed] , (37) 

Xo = Xo(Q,o) 

= 2 [p2(5~^)p„po +p^(5"^)p^_p^ - 2poP2(5"^)popJ . 

(38) 

The inversion of the 7x7 matrix is performed numeri- 
cally. 

After analytic continuation [iu^m — > w + iO+) the exci- 
tonic susceptibility yields 



xx(q,w) = 



(0)/ \ 

xx (q^w) 



xy°^(q,^) + i 



(39) 
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with 



and 



k 



kt 



Sri 



Pi 



2 



'e(O) 



k+q^ 



Po 



(40) 



(41) 



For the BI at T = the random phase approximation 

fi s 
result^ is recovered, — = U. 



To determine the EI phase we compute the static exci- 
tonic susceptibihty xx(q, 0). The direct band gap situa- 
tion gives the order vector of the EI phase as q = 0. Using 
Eq. ([TO]) the fluctuation matrix element Sp^p^ [Eq. PT|)] 
reduces to 



Pz 



(42) 



It is easy to see that the condition for the divergence of 
Xa'(0,0) equates to the gap equation ([28]) . 

The poles of Re Xx^ (Q: ^) [Eq- (|40)) ] give the continuum 
of electron- hole excitations, i.e., Wfc(q) = fi'k-i-qi — ^'kt- 
Excitonic pairing of electrons and holes is described by 
the pole of Rexx(q, w) [Eq. outside the electron- 

hole continuum,— i.e., by 



Pi 



(43) 



with 



where wc(q) 
given by 



<ojx{q) <uJci^) , (44) 
'^A;(q)|min- The cxcitou binding energy is 



-E|(q) = wc(q) -wx(q) 



(45) 



We want to emphasize that ujx , wc and E^ are explicitly 
q-dependent in contrast to Ref. |6|, where only excitons 
with q = are considered, and Ref. H, where the exciton 
binding energy is assumed to be independent of q. 



III. NUMERICAL RESULTS 
A. Instabilities against CDW and SOO 

To obtain results for the 3D EFKM we transform the 
k-summation into an energy integral using the tight- 
binding density of states (DOS) for a simple cubic lattice. 

From the charge and orbital susceptibility we derive 
information about the CDW and SOO formation, respec- 
tively. For asymmetric bands (|t-|-| \ti\) the charge and 
orbital susceptibility diverge at the same critical Ef, as 
shown in Fig. [21 implying \5^\ ^ \Si\. The analogy be- 
tween CDW and SOO vanishes if the bandwidths are 



equal, as can be seen for t^- = —1.0 in Fig. [2l In this 
case, the orbital susceptibility diverges contrary to the 
charge susceptibility, thus, a CDW will not develop and 
6^ = —Si- We conclude that the density inhomogeneity 
(5cr is largely affected by the bandwidth. 

Figure |3] shows Xo and Xc for = —0.8. The suscepti- 
bilities diverge at the same critical E^. With increasing 
strength of the Coulomb interaction the critical |£'-t-| for 
CDW (SOO) formation increases, because for a larger 
interaction the charge (orbital) order becomes more fa- 
vorable. Figure [3] clearly shows that the CDW and SOO 
region is confined close to the symmetric case E-^ = 0. 

For smah band splitting either the CDW (SOO) or 
the EI, separated by a first-order phase transition, can 
be realized, and one has to compare the free energies to 
identify the true ground state. Hence, to determine the 
SB ground-state phase diagram (analogous to the HF 
case shown in Fig. [1]) the generalization of the saddle- 
point equations to a bipartite lattice is inevitable, which 
is beyond the scope of this work. To investigate the EI in 
the following, we choose the band-structure parameters 
E^ = -2.4 and = -0.8, where a CDW (SOO) is not 
realized (sec Fig[3|). 



B. Instability against EI 

Figure |4] shows that the EI phase boundary in the 
weak-coupling as well as in the strong-coupling regime 
is reproduced by poles of the uniform static excitonic 
susceptibility, as demonstrated analytically in Sec. IIIDI 
To determine the region where free excitons can exist, 
we evaluate the condition for exciton formation (|43p sub- 
jected to the constraint The exciton binding en- 
ergy has to be positive. For numerical reasons we set 
the threshold to min(iJ^) = 10~^. For the 3D case we 
restrict ourselves to excitons with a center-of-mass mo- 
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FIG. 2: (Color online) Static charge and orbital susceptibility 
of the 3D EFKM for T = and (7 = 4 as a function of Ef. 



6 



1 - 



100- 







-100 



U = 2 

U = 4 

U = 6 



-1.5 -1 



-0.5 



FIG. 3: (Color online) Static charge and orbital susceptibility 
of the 3D EFKM for T = and = -0.8 as a function of 
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FIG. 4: (Color online) EI phase boundary (black solid line) of 
the 3D EFKM in the U-T plane. The red diamonds give the 
poles of the uniform static excitonic susceptibility for some 
fixed values of U. The blue circles give the critical Ux for 
exciton formation with center-of-mass momentum q = 0. The 
black dashed line separates the SM and the SC phase. Inside 
the EI phase we suggest that the black dashed line gives a 
crude estimate for the BCS-BEC crossover region (see text). 
The band-structure parameters are Ef — —2.4 and = —0.8. 



mentum q = 0. In Fig. |4] the critical Coulomb interac- 
tion strength for the exciton formation Ux{T) is shown 
by blue circles, where zero-momentum excitonic states 
exist for U > Ux- For low temperatures Ux{T) coin- 
cides with the EI phase boundary in the strong-coupling 
regime. This suggests that in this regime the EI is de- 
veloped by a BEC of preformed cxcitons. Within our 
analysis, for high temperatures Ux{T) deviates slightly 
from the critical Ug{T), determined from Eq. ((27)) . which 
separates the SM {U < Ug) and the SC {U > Ug), except 
for the point where Ug{T) coincides with the EI phase 
boundary. This deviation turns out to be a result of the 



required finite binding energy of the excitons. In a SM 
cxcitons with q = cannot exist. Here, the EI develops 
due to a BCS-type pairing of electrons and holes. We 
suggest that the BCS-BEC crossover region is estimated 
by extrapolating Ug{T) into the EI phase. 

To consider excitons with an arbitrary center-of-mass 
momentum, the q-resolved excitonic susceptibility is cal- 
culated for the 2D EFKM, in order to keep the numerical 
effort manageable. Again we choose the band-structure 



parameters E- 



-2.4 and t 



-0.8, for which the EI 



phase is stable in 2D )^^i^° see Fig. [T] To evaluate the 
SB parameters, the k-summation is transformed into an 
energy integral using the tight-binding DOS for a square 
lattice. For the computation of the excitonic susceptibil- 
ity the k-summation is explicitly performed. 
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FIG. 5: Phase diagram of the 2D EFKM for the band- 
structure parameters — —2.4 and — —0.8. The solid 
line shows the EI phase boundary and the dashed line sep- 
arates the SM and the SC phase. The exciton dispersion at 
the marked points a,b,c,d is shown in Fig. [HI and the exciton 
dispersion for e and f is shown in Fig. |9l 

The phase diagram is shown in Fig [5] Qualitatively 
there is no difference between the phase diagram of the 
2D and 3D EFKM. Hence, the following results for 2D 
should hold qualitatively for 3D, too. 

Figure [6] shows the static excitonic susceptibility for 
zero temperature. The formation of the EI phase is in- 
dicated by the divergence of xx(q, 0) at q = for the 
critical value Uei = 5.71. For U > Uei the static exci- 
tonic susceptibility remains finite for all q. 



C. Excitonic bound states 

Next we want to characterize the paraphase in the 
vicinity of the EI phase with a view to the formation of 
uncondensed excitons. Figure [7] shows the real part of the 
dynamic excitonic susceptibility outside the continuum 
for several values of tu. The divergences of Rcxx{<i,^) 
point out that the system is unstable against the for- 
mation of free cxcitons. With increasing excitation en- 
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FIG. 6: (Color online) Static excitonic susceptibility for T = 
as a function of q (2D). For all U we obtain — 1 and = 0. 



ergy w the exciton momentum q for the exciton forma- 
tion increases due to the direct band gap situation. The 



divergence of Kcx'xi'iT^) means 



ujc (q) , shown 



as the vertical solid lines in Fig. [T] For the q values 
where Re xx (q, ^) (upper panel) and Re Xx^ w) (lower 
panel) is not plotted in Fig. [7l the given uj is larger than 
wc'(q). Hence, these q values are irrelevant for the exci- 
ton formation for the considered excitation energy w. 

Having confirmed the existence of excitonic states, we 
now turn to the properties of these states. We denote 
the energy-momentum relation derived from Eq. (|43p as 
the dispersion of the exciton band. The smallest pole 
of Re Xx' (l' ^) defines the q-resolved lower boundary of 
the electron-hole excitation continuum, which we denote 
in the following as the continuum band. In Fig. [5] the 
exciton band ujx (q) and the continuum band ujc (q) are 
shown for characteristic points in the phase diagram (see 
Fig- El- In general, the continuum band is more sensitive 
to the temperature and Coulomb strength than the ex- 
citon band. Hence, the binding energy of the excitons is 
mainly affected by the continuum band when T or U is 
varied. Figure [5] suggests that the exciton dispersion has 
a cosine-like form, determined by the underlying lattice. 

In Fig. [Hlja), for any momentum a finite energy is 
needed to transfer one electron from the valence band 
to the conduction band, wc(q) > 0, which is character- 
istic for the underlying SC band structure. Likewise the 
creation of an exciton requires energy, ujx (q) > for all 
q. The comparison of Fig. ^h) with Fig. [SJa) unveils 
the influence of the Coulomb interaction strength. In 
Fig. [S^b), the exciton band is shifted to lower energies 
because the point (b) in the phase diagram (see Fig.[5|) is 
closer to the EI phase than (a). For the continuum bands 

w^^(q) < uj'^\q) and, therefore, the binding energy of 
the excitons in Fig. ^h) is smaller than in Fig. ^a) , since 
C/(b) < C^(a): i-e., the Coulomb attraction between elec- 
trons and holes in Fig.lSJb) is smaller than in Fig.[5Ja) as 
well, and the electrons and holes are more loosely bound. 
The exciton and continuum dispersion at the SC-SM 
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FIG. 7: (Color online) Real part of the dynamic excitonic sus- 
ceptibility as a function of q (2D) for U = 5.5 and T — 0.256 
(upper panel) and the corresponding ReXx'('3i'^) (lower 
panel). The densities are = 0.994 and = 0.006. The 
vertical solid lines show the lower boundary of the electron- 
hole excitation continuum. The black dotted line gives the 
critical value of Re^x^ foi" the exciton formation. 



transition are shown in Fig. [H^e) . The continuum band 
reaches w = for q = 0, due to the direct band gap 
situation. The excitonic band disappears for small, finite 
center-of-mass momenta, not only for q = 0. The reason 
is the required finite binding energy of the excitons, i.e., 
(q) should be separated from a;c(q). 

In Fig. IHJd) the SM band structure is characterized by 
the vanishing of positive ujc (q) at some finite value of q. 
In a SM band structure excitonic states exist only with 
finite center-of-mass momenta, in contrast to Figs, 
and^h). The existence of excitonic bound states in a 
SM is confirmed by a very recent EFKM study with the 
projector-based renormalization method»^ The compar- 
ison of Figs, mjb), M,c), andlHJd) demonstrates that the 
energetic position of the continuum band decreases with 
increasing temperature and, as a result, the exciton bind- 
ing energy is lowered. 

The qualitatively different feature of the exciton and 
the continuum band in the SC and in the SM phase sug- 
gests that the condensation process at the SC-EI transi- 
tion differs from the one at the SM-EI transition. 
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(a) T = 0.26 
U = 6.00 




(b) T = 0.26 
U = 5.50 



(c) T = 0.34 
U = 5.50 



(d) T = 0.40 
U = 5.50 



q /it 



q /it 



FIG. 8: Exciton band (solid line) in comparison with the 
continuum band (dashed line) at the points marked in Fig. (5] 

(a) in the SC phase {n^ = 0.996, = 0.004, p. = 0.527), 

(b) in the SC phase with a smaller band gap {rif = 0.988, 
nj. = 0.012, fl = 0.188), (c) at the SC-SM transition (nt = 
0.973, ni = 0.027, p. = 0.056), (d) in the SM phase (n^ = 
0.960, ni = 0.040, p = -0.013). The chemical potentials 
are measured relative to the valence band top. In all figures 
qy = (2D). 



(a) T = 0.20 
U = 5.50 




(b) T = 0.26 
U = 4.80 
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FIG. 9: Exciton band (solid line) and continuum band 
(dashed line) (a) at point (e) marked in Fig. [5] the SC-EI 
transition (n^ = 0.995, = 0.005, p, = 0.280), and (b) at 
point(D marked in Fig. [5] the SM-EI transition {n^ = 0.949, 
ni = 0.051, — —0.375). The chemical potentials are mea- 
sured relative to the valence band top. In all figures qy — 0. 



Figure inja) shows the exciton and the continuum dis- 
persion at the SC-EI transition. As already stated, 
^c{<i) > Vq reflects the underlying SC band structure. 
The minimum of the exciton band (at q = 0) reaches 
cj = 0, so that the creation of a zero-momentum exci- 
ton does not require energy. Physically, comparing only 



the energies the exciton band touches the valence band 
top. In our analysis, this is the point where the BEC of 
excitons sets in. 

On the contrary, the exciton dispersion at the SM-EI 
transition, shown in Fig. IHl^b), gives no hint for a con- 
densation process. In fact, it is a characteristic plot for 
a SM band structure: there are excitonic bound states, 
but only with a finite center-of-mass momentum, c.f. 
Fig. El^d). In the SM regime there is a large and well- 
defined Fermi surface and the electron-hole condensation 
at the SM-EI transition can be imagined as a BCS-type 
process. 



IV. SUMMARY 

In this paper we studied the extended Falicov-Kimball 
model within a S'0(2)-invariant slave-boson treatment 
taking Gaussian fluctuations into account. We computed 
the static charge and orbital susceptibility in order to in- 
vestigate the formation of a charge density wave and stag- 
gered orbital order. Analyzing the static excitonic sus- 
ceptibility, the instability towards an excitonic insulator 
(EI) phase agrees with the saddle-point phase diagram. 
By calculating the dynamic excitonic susceptibility, we 
analyzed the formation of excitons in the phase without 
long-range order. We found that finite-momentum ex- 
citons form in the vicinity of the EI phase, not only at 
the semiconducting (SC) side but also at the semimetal 
(SM) side. At the transition to the EI phase we ob- 
served qualitatively different features at the SC and the 
SM side. At the SC-EI transition the excitation energy 
of the excitonic state with zero center-of-mass momen- 
tum vanishes, leading to a Bose-Einstein condensation of 
these preformed excitons. At the SM side there are no 
excitonic bound states with q = 0. Here, the EI forms 
due to a BCS-type pairing of electrons and holes, and the 
occurring excitonic states of finite momentum are not of 
importance for the phase transition. 
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